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Using a flowchart representation of quantum optomechanical dynamics, we design coherent 
quantum-noise-cancellation schemes that can eliminate the back-action noise induced by radiation 
pressure at all frequencies and thus overcome the standard quantum limit of force sensing. The 
proposed schemes can be regarded as novel examples of coherent feedforward quantum control. 
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Noise introduces error to communication or sensing 
because it is unknown to the observer and cannot be 
distinguished from the desired signal. If the noise can 
be measured separately, one can physically or compu- 
tationally remove the noise. This is the principle of 
noise cancellation, a technique that has seen widespread 
commercial application, especially for acoustic-noise con- 
trol [H: noise-cancelling headphones, for example, work 
by recording the ambient noise and playing it back with 
opposite amplitude to interfere destructively with the 
noise reaching the inner ear. 

Here we consider the application of noise cancellation 
to quantum systems. We focus on optomechanical force 
sensors, in which quantum radiation-pressure fluctua- 
tions on a moving mirror within an optical cavity can 
introduce excess noise, called the back-action noise, in 
addition to the shot noise at the cavity output. This 
leads to the standard quantum limit (SQL) on detect- 
ing a classical force on the mirror Various methods 
of quantum-noise reduction have been proposed to over- 
come the SQL, including frequency-dependent squeezing 
of the input light , variational measurement [1, 0, 1| ; 
introducing a Kerr medium inside the optical cavity 9j, 
and the use of dual mechanical resonators 13, lli or an 
optical spring 12, 13] to modify the mechanical response 
function. Several of these proposals have been imple- 
mented 0, El, 13 1, though not in the quantum regime. All 
can be regarded as quantum-noise-cancellation (QNC) 
schemes, which utilize destructive interference to reduce 
or eliminate the effects of back action. With back action 
tamed, squeezing of the input light [H] can be used to 
improve sensitivity further. 

To facilitate understanding of the QNC concept and 
design of new QNC schemes, we introduce the use of 
flowcharts to depict the quantum dynamics. Using 
the flowcharts, we design a few novel QNC schemes 
that can eliminate back-action noise at all frequen- 
cies and should require less space than previous broad- 
band QNC schemes. Although the effects of quan- 
tum radiation-pressure noise are only now becoming de- 
tectable [15], the field of optomechanics has recently 
seen rapid progress [l6|. Back-action noise is expected 
to become a major issue in future optomechanical sen- 



Combined with quantum filtering and smooth- 
ing pj|, QNC has the potential to improve significantly 
the performance of future force sensors beyond conven- 
tion al q uantum limits. In the context of quantum con- 
trol [18[, QNC schemes can be regarded as examples of 
coherent feedforward quantum control, to be contrasted 
with measurement-based feedforward control [l9j and co- 
herent feedback control [13] techniques. 




FIG. 1: (Color online) Optomechanical force sensor. 

Consider now a Fabry-Perot cavity with a moving mir- 
ror, as depicted in Fig.[TJ The mirror is modeled optically 
as perfectly reflecting and mechanically as a harmonic os- 
cillator with position operator q{t), momentum operator 
p{t), mass m, and resonant frequency w„i, subject to a 
force f{t). The cavity is pumped with an input beam 
ylin(i), incident on the partially transmitting mirror on 
the left. We work in a rotating frame that removes the 
harmonic time dependence at the input beam's carrier 
frequency ujq. All phases are referenced to the input field, 
which has constant (real) mean amplitude {A-m) — Ain- 
The intra-cavity field decays at rate 7 due to coupling 
through the partially transmitting mirror to the output 
beam Aout(i)- The intra-cavity optical field a(t) is as- 
sumed to be resonant at the carrier frequency loq when 
the mirror is displaced to its equilibrium position under 
the mean radiation pressure, at which point the round- 
trip length is 2L. The mirror position q{t) is defined rel- 
ative to this equilibrium position. The force is estimated 
by continuous homodyne measurement of an appropriate 
quadrature of Aout, labeled 2 in the following and usu- 
ally the phase quadrature. The operators obey canonical 
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commutation relations: [q{t),p{t)] = ih, [a{t),a^{t)] = 1, 
and [A,^{t),Alit')]=5{t-t'). 

One can generalize this basic setup to more elaborate 
configurations, e.g., more complicated optical and me- 
chanical mode structures or detuned cavity excitation, 
which gives rise to an optical spring 12|, [l3|. A ma- 
jor difference from a two-armed interferometer is that 
squeezed light must be input in the same beam as the 
mean field that powers the system. Nonetheless, this 
minimal setup already captures the salient features of 
cavity optomechanical systems. 

The operator equations of motion are 
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where a = Ain\/2/^ is the mean cavity field ((a^a) = 
a^). Removing mean fields and defining amplitude and 
phase quadrature operators for the fluctuations that re- 
main, ^in = Ain + (Cl + i£,2)/V^, Aout = Ain + {Vl + 

ir]2)/V^, and a = a + (oi -I- ia2)/V2, one can linearize 
Eqs. ^ by neglecting quadratic terms and obtain a sys- 
tem of linear differential equations. 



dx 



= Fx + Gw , y = Hx + Jw 
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where n = y/2aijjQ/L is an optomechanical coupling 
strength. 

In control-system design, it is often illuminating to 
draw a block diagram to represent a system of differential 
equations [2lj . Here we use a simpler depiction, which we 
call a flowchart, as shown in Fig.[5Ja) for Eqs. © and 
Arrows point from a variable on the right-hand side of an 
equation of motion to a connected variable on the left- 
hand side. Since the system is linear, the output signals 
are sums of independent contributions from the inputs, so 
one can easily depict the flow of signal and noise from in- 
put to output. In Fig.[5][a), the dashed arrows depict con- 
tributions to the output quadrature 772 from the signal / 




FIG. 2: (Color online) (a) Flowchart representation of 
Eqs. ((2]) and (|3}- (b) Simplified fiowchart depicting only the 
input and output signals linked by transfer functions. 



and from the input amplitude and phase fluctuations, ^1 
and ^2- The ^2 contribution is commonly known as shot 
noise, while the ^1 contribution is the back-action noise 
due to the Kerr-like ponderomotive coupling of the cavity 
amplitude quadrature ai to the phase quadrature 02. 

The solution for the state variables contains a tran- 
sient solution, which decays exponentially and which can 
be ignored by pushing the initial time to —00. It is conve- 
nient to Fourier transform the remaining inhomogenous 
solution, writing y{Vl) = K{n)w{n), where the overall 
transfer matrix is 

K{n) = -Hm + F)-'G + J=(^ ° 

(4) 

with I being the identity matrix. The cavity and pon- 
deromotive transfer functions are 

(0\ - 7 + ^^ (c,\ - 27n^Vm 



(5) 



and the signal transfer function is Kf{^) = Kpni{^){l — 
ifj) / h,K,\f2Fi. A simplifled flowchart that connects outputs 
to inputs by transfer functions is shown in Fig. [2jb). 





FIG. 3: (Color online) (a) Frequency-dependent input squeez- 
ing, (b) Variational-measurement scheme. 

Unruh's proposal for frequency-dependent in- 
put squeezing [3|-|7[ and Vyatchanin and Matsko's 
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variational-measurement scheme 0, H [1| can be under- 
stood using flowcharts. First consider Unruh's proposal, 
flowcharted in Fig. [21[a). The input beam, with quadra- 
tures xi and X2, is processed through a passive system 
that rotates the quadratures by a frequency-dependent 
angle after which the beam is displaced by 

to produce the input field to the optical cavity. 
The transfer functions from xi a-nd X2 to the output 
phase quadrature 772 are K^^^xi — cos + i^cav sine' 
and JsT^^vo = - K„m sin 6 + Kc, 



^cav ^ 

L^2X2 — — ^ipm fji^r -"^cav cos 0. The rotation 
is chosen to make the contribution from xi vanish, 
i.e., i^,)2Xi ~ 0, by introducing destructively interfering 
paths from xi to 772, as shown in Fig.[3lja). This requires 
tan 0(17) = -Kp^{Q)/Kc^^{n). Kimble et al. devised 
a method for performing the dispersive rotation over a 
wide bandwidth by filtering the input beam through two 
detuned cavities. If the input to the filtering cavities 
is vacuum, the rotation has no effect on sensitivity; 
since X2 is now entirely responsible for the output noise, 
however, it can be squeezed to improve the sensitivity 
to the shot-noise limit and beyond. 

The variational-measurement scheme can be under- 
stood using the flowchart in Fig. |31[b) . This scheme dis- 
persively rotates the output light and measures the fi- 
nal output quadrature C2, which is a combination of the 
original output quadratures ryi and 772- This introduces 
an "anti-noise" path from ^1 to (,2 via 771, which can be 
used to cancel the original back-action noise path via 
r]2- Given the similarity of the flowcharts in Fig. [3l the 
required dispersive rotation is the same as that for Un- 
ruh's proposal and can again be implemented using two 
detuned cavities. With back-action noise eliminated, the 
sensitivity is limited by shot noise and can be improved 
further by squeezing the input quadrature ^2 @]- 

While the aforementioned schemes can be regarded as 
examples of coherent QNC, it is clear from Fig. [2] that a 
more direct way of cancelling the back-action noise is to 
introduce an anti- noise path from ai to 02 , as illustrated 
in Fig. m This calls for coherent processing of the intra- 
cavity field, using parametric interactions to undo the 
ponderomotive squeezing. The use of a Kerr medium Q , 
the dual-mechanical- resonator scheme 13, III, and even 
an optical spring 12 , 13] can also be thought of as intra- 
cavity QNC schemes, but they cannot eliminate the back- 
action noise at all frequencies. To achieve broadband 
QNC, we introduce two auxiliary state variables, q' and 
p', which play the role of q and p in the anti-noise path, 
as shown in Fig. HJc). Adjustable parameters 5, ^, and 
v, all with units of frequency, characterize the couplings 
in the anti-noise path. The anti-noise transfer func- 
tion, indicated in Fig. ^h) by the green dash-dot ar- 
row, is 2jfig'^ /{—fiv — r2^)(7 — ifi)^ — — i^pm- We need 
fiiy — —uj'^ and /ig^ = —hK^/m, thus requiring ^ to be 
negative. This cannot be implemented by ponderomotive 
coupling to another mechanical oscillator unless its mass 
is negative. 
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FIG. 4: (Color online) (a) Back-action noise cancellation 
achieved by introducing an anti-noise path (green dash-dot 
arrow) from the intra-cavity amplitude quadrature ai to the 
phase quadrature 02. (b) Simplified fiowchart. (c) Detailed 
flowchart of the ponderomotive coupling and the intra-cavity 
matched squeezing, (d) An implementation of the matched 
squeezing scheme. 



The matched squeezing can nonetheless be achieved 
by coupling the intra-cavity field to an auxiliary field b = 
{q' + ip') I \pl. The required parameters are v = = 
u!,n and g = Ky^h/mLu„i — au!Qy^2h/'muJm/L, and the 
equations for a and b become 



da 
db 



— — —7a - 



—qa + —{b + W) 



— = iLOnib + —{a + flt) - iga 



(6) 



The auxiliary field should be inside another optical cav- 
ity, with resonant frequency wq — w™; it plays the role 
of a negative- energy mode in the anti-noise path. The 
constant driving term —iga, which removes the mean 
field from 6, can be eliminated by redefining the auxil- 
iary mode as 5' = 6 — ga/ujm and displacing the input 
field Ain, resulting in the following equations of motion: 



^ga+|(6' + 6't) + V2^A;„ 



da 

dl^-^" 

db' . , , ig , |. 
— = lUjnb + -j(a + a' ) , 



AL = Au 



(7) 



The coupling between a and b' can then be realized using 
a beamsplitter (BS) and an optical parametric amplifier 
(OPA), as schematically shown in Fig. Sfd). With the 
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back-action noise removed, only the shot noise from ^2 
remains, and one can improve sensitivity further by in- 
creasing the optical power or by squeezing ^2- 

In terms of the coupling constant g, the ponderomotive 
transfer function has the form 



(8) 



so |i4rpm(fi)|/|-ftrcav(f^)| ~ Q^ll'^m for frequencies away 
from the mechanical resonance; back action is thus im- 
portant when g > ^/^^uJ^■ The required single-pass idler 
gain of the OPA is thus G ~ {gL/cf > 7a;,„ (L/c)^ = Gq- 
Go ranges from about {2-k x 10 Hz x 5 km/c)^ ~ 10^^ 
for gravitational-wave detectors to {2tt x 50 MHz x 
100 pLTa/cY — 10~^ for microscale systems [l^, so the 
required G should be easily achievable with current OPA 
technology. 

There is another potential problem: the number of 
photons in the auxiliary cavity in Fig. Hfd) is {b'^b') — 
g^o? jujl^^ which becomes significantly higher than that 
in the primary cavity when a;„i is small, making the intra- 
cavity scheme problematic for applications that use high 
circulating power. This problem can be alleviated by 
using matched squeezing to modify the input or out- 
put optics. This requires the use of another double- 
cavity setup, similar to that in Fig. Hl^d), but without 
the moving mirror, to pre-squeeze the input light go- 
ing into the sensor cavity or post-squeeze the output 
light, as shown in Figs.[5{a) and (b). The flowcharts for 
these schemes, shown in Figs.[5jc) and (d), demonstrate 
broadband cancellation of back-action noise much like 
the intra-cavity scheme. If one is interested only in low 
frequencies $7 <C 7, the input/output matched squeezing 
can be implemented using smaller cavities with a larger 
decay rate 7 and a larger coupling constant with 5^/7 
held constant. 

We have assumed an ideal system, neglecting intrinsic 
mechanical and optical losses, to illustrate the essential 
features of QNC. Such assumptions are unrealistic in 
practice. To assess QNC schemes, it will be important 
to take into account the fluctuations associated with dis- 
sipation, such as the thermal noise associated with me- 
chanical damping. The design and performance of QNC 
in the presence of realistic dissipation and noise deserve 
further investigation. 
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